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Abstract 

Using Komar's definition, we give expressions for the mass and 
angular momentum of a rotating acoustic black hole. We show that 
the mass and angular momentum so defined, obey the equilibrium 
version of the first law of Black Hole thermodynamics. We also show 
that when a phonon passes by a vortex with a sink, its trajectory 
is bent. The angle of bending of the sound wave to leading order is 
quadratic in A/cb and B/cb, where b is the impact parameter and A 
and B are the parameters in the velocity of the fluid flow. The time 
delay in the propagation of sound wave which to first order depends 
only on B/c 2 and is independent of A. 

1 Introduction 

In 1981, Unruhp!J|2] came up with the idea that the phenomenon of Hawking 
radiation could be observed in an inviscid fluid which has barotropic equation 
of state at sub-Planckian energy scales. The idea for this model arose from 
a very simple question as to how a sound wave propagates in an inhomoge- 
neous fluid. He showed that if the fluid is barotropic and inviscid, and the 
flow of the fluid is irrotational, the equation of motion that the fluctuation 
of the velocity potential obeys is identical to that of a minimally coupled 
massless scalar field in curved space-time Lorentzian geometry. Since the 



fluid has barotropic equation of state, pressure force does not generate vor- 
ticity. Thus an irrotational flow will remain ir rotational. The viscosity of the 
fluid is ignored in order to maintain the Lorentzian character of the acoustic 
geometry. The fluid-flow will contain an analogue of a black hole if there is 
region where it flows with supersonic speed. The boundary of this acoustic 
black hole is a surface where the radial velocity of the fluid exceeds the local 
speed of sound. It has also been shown by Unruh that an acoustic black 
hole emits thermal radiation at a temperature given by the surface gravity 
at the sonic horizon. In his demonstration he has chosen a perfect sink flow 
(v = —A/r r) of the fluid which corresponds to non-rotating black hole at 
r = A/c in the acoustic geometry. Later this idea was extended by Matt 
Visser to a more general situation where he considered a fluid motion which 
has both radial and azimuthal components. This fluid flow corresponds to 
rotating black hole at the same position in the acoustic geometry. This cor- 
respondence has been borne out by the fact that such an acoustic rotational 
black hole exhibits an analogue of superradiance as shown by Basak and 
Majumdar [HUE]- As there is no analogue of the Einstein equation that the 
fluid under consideration obeys 1 , it can mimic only the kinematical aspects 
of gravity, not its dynamics. 

For asymptotically flat solutions of Einstein's equation one can define the 
Komar charges corresponding to each isometry of the spacetime. In the case 
of spacetime around a Kerr black hole Komar charges corresponding to time 
translational and rotational Killing vectors are called the mass and angular 
momentum of the black hole and they are related by the equilibrium version 
of the first law of black hole thermodynamics. The basic criteria, for defining 
Komar charges and also for those charges to obey equilibrium version of the 
first law of black hole thermodynamics, have been discussed in section(3). 
Interestingly acoustic geometry can also satisfy all those criteria. Hence we 
can define conserved charges for acoustic black hole . Since the equilibrium 
version of the first law of black hole thermodynamics does not depend on Ein- 
stein's equation, Komar charges for the acoustic black hole will also satisfy 
this law. 

One of the remarkable predictions of general relativity is that light tra- 
jectories are bent in passing by a massive object because the spacetime in 
the neighborhood of that massive object is curved. The bending of light in 
the neighborhood of a massive object can be described to sufficient accu- 

lr The fluid motion is governed by the Euler equation and the continuity equation 
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racy using kinematical considerations alone since the effect of the light on 
the geometry is negligible. Moreover the source of curvature is not directly 
relevant to the discussion; the light bending effect is a direct consequence of 
the curved geometry which is locally Lorentzian. Since acoustic geometry is 
an example of Lorentzian geometry, it is expected that this phenomenon will 
also occur in this geometry if we consider a curved geometry arising from a 
suitable choice of velocity of the fluid flow. 

Based on the idea that sound waves experience an effective gravitational 
field when passing through an inhomogeneous vorticity-free flow of a barotropic 
inviscid fluid, Fischer and Visser jH] showed that when a sound wave passes 
by an irrotational vortex, its trajectory is bent. This is analogous to the 
light bending effect in a physical gravitational field. In this paper we have 
demonstrated, using differential geometric techniques which we normally use 
in General Relativity, the bending of a sound wave for a spiral flow ( vortex 
with sink ) of the fluid. We have assumed that the background density of the 
fluid is constant, which implies that background pressure and the speed of 
sound are also constant and we have treated the sound wave as a ray. So our 
analysis is valid when is very less compared to the typical scale of varia- 
tion of the acoustic geometry, but it should be greater than the interatomic 
distance because we have considered the fluid as a continuous medium. 

The plan of the paper is as follows. We first briefly discuss in section (J2J), 
the geometry of an acoustic black hole. Using the geometric properties of 
acoustic space-time, we define the mass and angular momentum of an acous- 
tic black hole in section (jSJ). In section (0J) we discuss the eikonal approxi- 
mation which we will use to study the propagation of a sound wave in this 
geometry. In section (jSj) we present an analysis of null geodesies in acoustic 
geometry. This is followed by a demonstration of the bending of a phonon in 
acoustic geometry and the corresponding time delay in section (JHJ) and (J7j) re- 
spectively. Experimental aspects have been discussed in section (jSJ). Finally 
we conclude with a summary of the results obtained in this paper. 

2 Acoustic Geometry 

We consider the draining bathtub type of fluid motion with a line of sink at 
the center, which is basically a (2+1) dimensional flow with a sink at the 
origin. This leads to the velocity profile [7], [E] and [5] (in cylindrical polar 
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coordinates) , 



v = r H 



A . B 



(1) 



where A and i? are real and positive . 

A two surface, at r = — in this flow, where the fluid velocity is everywhere 
inward pointing and the radial component of the fluid velocity exceeds the 
local sound velocity everywhere, behaves as an outer trapped surface in this 
acoustic geometry. This surface can be identified with the future event hori- 
zon of the black hole. Since the fluid velocity is always inward pointing, 
the linearized fluctuations originated in the region bounded by the sonic hori- 
zon cannot cross this boundary . 

A linearized fluctuation in the dynamical quantities (density and pressure 
of the fluid, and velocity potential of the fluid motion) leads to the fluctua- 
tions of velocity potential and the equation of motion it satisfies is identical 
to that of a minimally coupled massless scalar field , 



If we assume that the background density of the fluid is constant, it auto- 
matically implies that background pressure and the local speed of sound are 
also constant. Thus we can ignore the position independent pre-factor in the 
metric because it will not effect the equation of motion of fluctuations of the 
velocity velocity potential. As for the Kerr black hole in general relativity, 
the radius of the boundary of ergosphere of an acoustic black hole is given 
by vanishing of goo> i- e > r ergo — VA 2 + B 2 /c. 



3 Mass and Angular Momentum of Acoustic 
Black Hole 



In this section I will discuss the definitions of certain conserved quantities 
in the case of a acoustic black hole which can be thought of as its mass 





where is acoustic metric and it is given by, 




2B dcp dt + dr 2 + r 2 d<p 2 (3) 



4 



and angular momentum. This is identical to Komar's definition of mass and 
angular momentum using isometries of spacetime. In his work, almost 40 
years ago, Komar|3] showed that for every isometry (£) of space time there 
exists a charge conserved(Q^) on spatial hyper-surfaces(S) which is given by, 

Qt = j[ X (4) 

This conserved charge is arbitrary up to a constant factor, which can be 
fixed using a known solution to the Einstein equation. For a Kerr black hole 
the conserved charges corresponding to the time translational and rotational 
Killing vectors are the mass(M) and angular momentum( J) of the black hole 
and they are given by, 

Mh = -t-^ / *dk ; J H = - J— / *dm (5) 

87T G jh 16n G Jh 

where k and m are time translational and rotational Killing vectors respec- 
tively. The signs in equation© reflect the signature of the spacetime geom- 
etry. 

Using this result Smarr^U] showed that if a stationary axisymmetric space- 
time contains a black hole then the mass and angular momentum of the black 
hole with respect to an observer at the horizon and one at asymptotic infinity 
are related as, 

m = m "-^gL' r ^ (6 > 
J = J * + ak L' R{m) (7) 

where *R{k) is Hodge dual of R{k){= R^^dx^). M and J are mass and 
angular momentum of the black hole with respect to a stationary observer at 
infinity, and Mh and Jh are the mass and angular momentum of the black 
hole with respect to an observer at horizon. Since for vacuum solutions of the 
Einstein's equation, R(k) and R(m) vanish, the mass and angular momentum 
of a black hole with respect to observers at the horizon and at asymptotic 
infinity are the same. If the spacetime contains a Killing horizon such that 
the normal to the horizon is a linear combination of time translational Killing 
vector and rotational Killing vector as, 

i = k + tt H m (8) 
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where f2# is the angular velocity of the horizon, then the mass at the horizon 
is given bv[TT].[T2]. 

M H = + 2 Q H Jh (9) 

If the black hole has electric charge then the above equation is modified to, 

M = ^r4?- + 2Q h Jh + V H Q H + .... (10) 

where Vh is the electric potential at the horizon and Qh is the charge of 
the black hole. If we restore all the occurrences of c, equation© is modified 
to the following form, 

M H = ^^ + 2 -^ (11) 
4 7r G Cr 

Even though in order to derive the equation (|T0|) for charged black hole, one 
has to use Einstein's equation with energy-momentum tensor of electromag- 
netic fields, equation (fTT|) is insensitive to the Einstein's equation. This equa- 
tion entirely depends upon the geometric properties of the spacetime and 
does not depend upon how geometry of the spacetime is changing with time. 
The above equation will hold if a curved spacetime has the following prop- 
erties, 

• Geometry of the spacetime is pseudo-Riemannian with signature (- 
,+,+,+) 

• Spacetime is asymptotically flat which implies that with increasing r , 
conserved current, = —R^ u ^ u decreases faster than the increase of 
area of the constant radius surface. 

• Killing vectors of the spacetime geometry are time-translational and 
rotational. 

• A Killing horizon is present with outward normal I = k + Qh m 

Acoustic geometry corresponding to the metric(JH} is similar to the geome- 
try of a curved spacetime containing a rotating black hole as far as kinematics 
are concerned. There are two isometries in this geometry and the correspond- 
ing Killing vectors are = 5f and m fI = 5^. The normal to the horizon 
is Z M = k^ + fi#m M . Since it is a linear combination of two Killing vectors 
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and fin = is constant, this vector is also a Killing vector of the acoustic 
metric. Therefore the event horizon of the acoustic black hole is a Killing 
horizon. So we can define the surface gravity of that Killing horizon as, 

PV^ = ^fl v (12) 

where Kh is the surface gravity of the horizon which is equal to c 3 /A. Hence 
this geometry is not only Lorentzian, but also satisfies all the four properties 
mentioned above. Therefore we can define the mass and angular momentum 
of the rotating acoustic black hole as the conserved charges corresponding to 
the time translational and rotational Killing vectors of the acoustic geometry. 
According to this, the conserved charges at the horizon corresponding to the 
time translational Killing vector is given by, 

Qk = ~ I *dk 
8 7r Jh 



'-[ d lx Ke^ a dx a 

(13) 
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(A 2 + B 2 ) c 



2 A 2 

and the conserved charge at the horizon corresponding to rotational Killing 
vector is given by, 

c 3 



Qm = —r- / *dm 
lb 7r Jh 



c 3 



16 n Jh 



Be 2 



d„m„ e dx a 



(14) 



vfe = — — 1 - 2 (15) 



4 

These two quantities satisfy the following relation , 

4 7T 6 



Here L H {= 2tt^) is the length of the horizon. 

This equation is identical to the relation between mass and angular momen- 
tum of physical rotating black hole in G = 1 unit. 
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The most interesting feature of these results is that the angular momentum 
is proportional to the velocity parameter B. If one switches off the rotational 
motion of the fluid, B will become zero and hence the angular momentum of 
the black hole disappears. This is expected because in the absence of B, the 
horizon is no longer rotating. If the vortex motion is quantized, i.e, if B takes 
values proportional to an integer, then Q m and Qk take discrete values. From 
the expression of Q m and Qk it is clear that Q m is proportional to integer and 
Qk has a minimum value c 2 /2 . So we can say that this acoustic black hole 
is analogous to a physical black hole whose mass and angular momentum 
are quantized. Since there is no analog of Einstein's equation in acoustic 
analog models of gravity, the relation between the variation of Qk and Q m 
is not identical to that of physical process version of first law of black hole 
thermodynamics. 

4 Eikonal Approximation 

We now analyze the geometric properties of acoustic space-time in terms of 
the null geodesies. We will consider the propagation of a sound wave in this 
geometry as a ray. This is valid in the eikonal approximation [TT].|13j which 
we describe briefly. 

The eikonal approximation is described by the two conditions, 

• Reduced wavelength (^) of the sound wave is very small compared to 
the typical radius of curvature of the acoustic geometry. 

• Reduced wavelength (^) of the sound wave is very small compared to 
typical length over which the amplitude, polarization and wavelength 
of the wave change. 

Under these conditions we can split the scalar field \& into a rapidly changing 
real phase and a slowly changing complex amplitude, 



Here e is a very small parameter. If R is the typical radius of curvature of the 
acoustic geometry and L is the typical length over which the amplitude and 




(16) 
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the polarization of the wave remains unchanged, then the order of magnitude 
of e is I Min(R, L) . The first term ( a ) in the amplitude is wave length 
independent. It is the dominant term in the expression while the rest of the 
terms are wave length dependent and also very small compared to the first 
term. Mathematically, the scalar field ( \1/ ) corresponds to a massless particle 
in the acoustic geometry when its frequency tends to infinityi.e, (\/2n — > 
). Practically all the particles have finite de Broglie wave length. This 
causes the wavelength dependent corrections to the amplitude of the scalar 
field. These correction terms disappear in the \/2tt — > limit. 

Substituting this in equation® and then collecting terms of order and 
- we get, 

V"0 = (17) 

and, 

V0 V M a = _~(V M V*0) a (18) 

Since sound rays are defined as the curves normal to the constant phase 8(x) 
surface and the wave vector = V M is normal to these surfaces, the tra- 
jectory of a sound ray is described by integral curves of wave vector. 
According to this definition of wave vector (k), 

k^ = (19) 

i.e, trajectories of sound ray are null curves and 

P V M a = -^(V M fc") a (20) 

Since the wave vector is the gradient of a scalar function, equation (fTUj) 
implies that the integral curves of the wave vector are also null geodesies i.e, 

k»V^k v = (21) 

Equation (j2(J|) leads to two important results, 
(1) If we define polarization of sound wave as, 



/aa 

then, 



/ = -4= (22) 
\Jaa 

/7=1 (23) 



9 



and, 



WVJ = 0, 



(24) 



So the polarization f is parallel transported along the integral curve of k^. 
(2) We can define a conserved vector j M = aak^ satisfying 

= (25) 

The corresponding conserved charge on spacelike hypersurfaces of acoustic 
geometry is, 

Q = jf f rf 3 S M (26) 

Physically this conservation corresponds to conservation of sound rays ( 
phonon number). Therefore in the eikonal approximation, we can treat the 
sound wave as a ray (phonon) which is a particle of zero rest mass with four 
momentum p = fik moving along null geodesies and parallel transporting the 
polarization f. Hence we will consider a phonon instead of sound wave in the 
analysis of null geodesies in acoustic geometry. 

The curvature (Ricci Scalar) of the acoustic geometry is, R = 2 r 2 ergo j r 4 . 
Our analysis is valid if ^- <C r 2 / 'y/2r ergo . Physically, when the wavelength of 
the phonon is less than the scale of variation of the metric, the effect of the 
curvature on the phonon's trajectory is significant. 



5 Null Geodesies in Acoustic Geometry 

As mentioned earlier there are two isometries (£a) of the acoustic geome- 
try. Corresponding to these we have the quantities conserved along the null 
geodesies given by, 

K A = -(UT 9,u ^ (27) 

Here A is the afline parameter of the geodesies and A takes values 1 and 
2. Now using the components of the acoustic metric in the above equation 
we get the following equations, 
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In order to find the expression for ^ we have made the following substitutions 
in equations, 

dt = dT+ r2 ^/_ A2 dr (30) 



A B 

d<p = dx+ , 2 2 t^t dr (31) 

r(r z c — A z ) 

Due to these transformations, equations ([28 )1 and (|29 p are modified to follow- 
ing equations, 



K2 = B dx~ r Tx (33) 



Solving equations (J52j) and (JBHj) we get, 

K x r 2 +K 2 B 



dX (r 2 c 2 - A 2 ) (34) 



rfx = Kj B r 2 + K 2 (A 2 + B 2 - r 2 c 2 ) 

dX ~ r\r 2 c 2 - A 2 ) { } 

Since A is the affine parameter of the null geodesies in acoustic geometry, 
the trajectory of a phonon can be determined by the following equation, 

dx^ dx v , , 



Using the components of the acoustic metric in the above equation, we get, 
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Substituting 4? and in equation (j3T|) we get, 

(dr\ 2 = K! 1K X K 2 B K\( A 2 + B 2 \ 

\d\) c 2 r 2 c 2 r 2 ^ r 2 c 2 J K } 

Rearranging the above equation we get, 

1 fd r Y ~\ t t 1 i^ n 2 



where, 



This equation describes the radial motion of a sound ray in acoustic geometry 
and it is identical to the one dimensional non-relativistic motion of a particle 
of unit mass and energy ^-^ i n a potential V(r) |TT]. 
The extrema of this potential, where dV ^ = 0, are at, 



= 2*,(A» + g) 

\ K 2 c 2 -2B K x y } 

There is a difference in the nature of the effective potential V(r) for < 
as compared with for K 2 . If K 2 < is negative then this potential has 
maxima at, 



2 K 2 (A 2 + B 2 ) 

for all negative values of K 2 because the 2nd derivative of this potential with 
respect to r is negative at r = tm for all negative values of K 2 . But when 
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K 2 > , V(r) does not have any extrema (i.e 



dV(r) 
' dr 



, does not have any 



real root in r) if K 2 c < 2 B K\. Therefore when K 2 is positive, for the 
existence of a maximum in V(r), we must have K 2 c 2 > 2 B K\. If a sound 
wave moves towards the center of the vortex such that V(tm) > it will 

be reflected at the turning point where 4v = . Considering the positive 
solutions of r of equation (jHHJ), we get position of the largest turning point of 
sound ray at, 



K 2 ( K 2 c 2 - 2 B Kj] 
2K, 2 



1 + 



_ 4 Ki 2 (A 2 + B 2 ) 
T ~ (K 2 c 2 - 2 B Ki) 2 



(43) 



If b(= \K 2 \ c/Ki) is the magnitude of the impact parameter and (A, B) <C cb 

_A_ 

c b 



then expanding in series (up to 2nd order) in the small quantities -A and 



A we get, 
r T = b 



1 - 



1 A 2 



2 c 2 6 2 



jj 2 
c 2 6 2 



+ O 




(44) 



This shows that, when K 2 is negative (i.e angular momentum of the sound 
ray is positive), tt is greater than impact parameter, b . Therefore in order 
to make the distance of closest approach shorter than the impact parameter 
we have to consider positive values of K 2 . Since 4| takes both positive 
and negative value, to find the bending angle we have to choose the sign 
of The choice is made on the basis of fact that when A = B = 0, we 
should obtain the Newtonian result, i.e, duration of propagation sound wave 
through the fluid should be positive. 



6 Bending of Phonon Trajectory 

In this section we will calculate the angle of bending of a phonon ray, ex- 
panding around A = and B = 0, as it passes by irrotational vortex. Before 
calculating the angle of bending one important thing has to be noticed. Since 
^ is not symmetric about the turning point even if ^ flips sign at that point, 
the change in the polar angle when the phonon travels towards the vortex is 
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not same as that when phonon travels away from the vortex. Therefore the 
total change of the polar angle throughout its path is given by, 



a , f rT dd) , r°° d<p 

A0 = / -f dr + / -f 

Joo dr Jr T dr 

= A0_ + A0+ 



d<f) 



dr 



(45) 



Expanding A0_ around A = B = we get, 



A 2 ( d2 /a, A B 2 ( d 2 , x , A 



+ A B 



d 2 



\dAdB 
-sgn( J fT 2 ) 



A=_B=0 
12 , D 2 



7T 3 7T (A 2 + £? 2 ) 

2 + 8 6 2 c 2 



- °m •(§) } 

Doing the same for A0 + we get, 

A0 + = (A0 + )a=b=o + A (^(A0+; 



1 A 5 



2 6 2 c 2 



A =J B=0 



(46) 



A=B=0 



A2 ( d \ AA A fi 2 / d 2 . A 



A=B=0 
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sgn(K 2 



7T 3_7T (A 2 + B 2 ) 

2 + ~8 



b 2 c 2 



+ 



1 A B 



2 b 2 c 2 



Hence the total change of polar angle of a phonon trajectory is, 



A</> = -sgn(iT 2 ) 



7T 



(47) 



(48) 



The first term corresponds to the case when fluid velocity v = 0. This 
matches with the Newtonian result. We have considered terms to leading 
order in A and B, which in this case is the 2nd order, but the terms corre- 
sponding to AB for A0_ and A0 + cancel each other and hence do not appear 
in the expression. Physically when the phonon approaches the turning point 
its motion is aided by the radial fluid motion, while it is hindered after cross- 
ing the turning point. The actual angle of bending is given by, 




This is consistent with Fischer and Visser'sffj] result for an irrotational vortex 
without a sink. Here the integrals are slightly tricky because the radius of 
the turning point depends both on A and B. So when we differentiate 
with respect to A and B , we take this into account. We have done the 
integration from ry to oo. This is a good approximation because the radius 
of the boundary of the vortex is very large compared to the turning point. 



7 Time delay and advance 

Unlike in the case of the bending of a phonon, in order to calculate the time 
delay and time advance [TH Hoi fTSl [T71 ITH1 ITU] in phonon propagation, we will 
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compare null geodesies of same turning point for the fluid with motion and 
the fluid at rest. Using equation we express K 2 in terms of r T , 



Substituting K 2 in ^ and 4| we get, time interval in which a phonon passes 
through the vortex, 

f r T dt , f r v dt , 

At = — dr + / — dr 

Jr v dr Jr T dr 

= At_ + At + (51) 
where ry is the radius of the vortex. 

Expanding At- around A = B = for fixed Tt we get, 

Ai - = {At - u = B =° + A 0->L M + B O^L. 

' A\ 2 /B\ 2 ) 



-yjr' v - 4 + 8gn(A" 2 ) ^tan" 1 h/(-V - 1 



Cr If Vry/ 

- HS) +0 {(l) 2 -(l) 2 } 

Doing the same for At + we get, 

Ai + = (Ai+) — + * (> +) L.„ + B (> t+) L. 

A\ 2 fB\ 2 ) 



- 4 + sgn(A- 2 ) | tan" 1 h/^)'-! 
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In contrast to the angle of bending, the time delay to leading order is linear 
in A and B. If we assume that (A/c 2 t, B /c 2 t) <C 1, where r = 1 



,we can neglect the higher order terms. Subtracting the zeroth order term 
from At we get the time delay up to first order in A and B to be, 



St 



At - At 



A=B=0 



IB 

= sgn(K 2 ) — tan -1 




+ o 



The only difference is in the limits of integration. This is because of the 
fact that the zeroth order term blows up if we put infinity in place of the 
radius of the vortex. In fact, the radius of the vortex, in actual experiment 
is large but finite. So a phonon will pass through the vortex in a finite time 
interval. The effect of the radial flow of the fluid on the propagation of 
the phonon is different during the two stages of motion - before and after 
having crossed the turning point. When the phonon approaches the turning 
point its motion is aided by the radial fluid motion, while it is hindered after 
crossing the turning point. The effect of the rotational velocity of the fluid 
is the same throughout the motion. Depending on the sign of K 2 and the 
angular momentum of the fluid, the sign of the time "delay" is decided by 
the combined effect of these. The terms that are first order in A/c 2 of At + 
and At_ cancel each other. This is because of the fact that when the radial 
coordinate is decreasing the radial flow of the fluid favors its motion and in 
the opposite case it opposes its motion. 



8 Experimental Outlook 

Unlike Hawking radiation and Superradiance for gravitational black-holes, 
the phenomena of bending of light and the corresponding time delay, as it 
passes by a massive object, have been experimentally observed. Hence exper- 
imental confirmation of the bending of the trajectory of a phonon will give 
further credence to our view that the phenomena of bending of the trajectory 
of a massless particle (and the corresponding time delay) in curved geome- 
tries, are features that are generic to Lorentzian geometry and independent 
of the dynamics of the system under consideration. In fact it would liberate 
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these phenomena from the almost stereotypical association with gravitating 
systems alone. For a given de Broglie wavelength of a phonon, there is a 
restriction on its impact parameter. If the impact parameter of a phonon is 
smaller than its de Broglie wavelength then the phonon behaves quantum me- 
chanically and hence both the concepts of trajectory and impact parameter 
are meaningless. On the other hand, classical analysis of phonon propagation 
can be done, when the de Broglie wavelength is very very small compared 
to the impact parameter. So, if we consider the semiclassical scattering of 
phonon in acoustic geometry, the minimum value of the impact parameter 
should be equal to the de Broglie wavelength of a phonon. Therefore the 
maximum value of the angle of bending (up to 2nd order) of the phonon 
trajectory and the corresponding time delay (up to 1st order) are given by, 



3 7T 



ergo 



A 



ph 



(55) 



St 



2B 



tan 1 




where X„h is the wave-length of phonon. 



(56) 



Since vortex motion with a sink has not yet been realized in the labora- 
tory, it is not possible to give a correct estimate of the angle of bending and 
the corresponding time delay, but it is certainly possible if we consider per- 
fect vortex flow(A=0). For a perfect vortex flow equations (jB3|) and (JEBj) are 
modified to, 

Q „ /„ \ 2 

(57) 



St 



2B 



tan 



3 7T / 


' ergo 


4 


f 








-l 




( ry 









(58) 



where r e 



_ B 
1 ergo — c ■ 

Instead of a single phonon if we consider two phonons propagating in 
the same direction but on opposite sides of the vortex, they will be seen to 
converge to a point. Thus a vortex in an inviscid, barotropic fluid acts like a 
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convergent lens and the corresponding focal length(up to 2nd order) is given 
by, 

2 A 3 

fmin = 7, — n (59) 

3?r r 2 ergo 

The circulation for a perfect vortex flow (v = B/r <fi) is defined as, 

T = j>v.dx = 2irB (60) 

If we consider a singly quantized vortex(i? = h/m), then the circulation of 
the vortex is,r = — . 

m 

As discussed in [§] there are two cases where we see perfect vortex motion, 

1. Superfluid Helium, 4 ife, 

2. Bose-Einstein condensation in atomic gases. 

(a) Bose-Einstein condensation of 23 Na, 

(b) Bose-Einstein condensation of 87 Rb, 

(1) Superfluid Helium, 4 He : In case of superfluid helium |2U] . |22|.|23]. 
[24J the energy spectrum of the quasiparticle, as shown by Feynman, is, 

» 2 * 2 , , 

<<*> = 2^ (61) 

where S(k) is the structure factor of the liquid which can be determined from 
neutron scattering. For small k, e(k) is linear and for k = 1.92 A , S(k) has a 
maximum . This corresponds to a minimum in e(k). The excitations around 
that minimum are called rotons and those corresponding to the linear part 
of the spectrum are called phonons. In case of superfluid helium the roton 
minimum occurs at wave length A = 3.27A and the spectrum is phononic in 
the wavelength range roughly three times of it. At temperature very close 
to 0°K , the speed of sound in helium is 240 m/s. So for a singly quantized 
vortex in liquid Helium, 

• Wavelength of phonon (X p h) = 9.8lA. 

• Radius of the ergo sphere, (r ergo ) = 0.655A. 

• Maximum angle of bending of phonon, 5<p = 0.6° 
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Minimum focal length,/ 



467A 



Maximum time delay is, St = 0.857ps. 



(2) Bose — Einstien Condensate : In case of a dilute Bose-Einstein con- 
densate |2Sl EH EHl Em EI3 E2| the energy spectrum of the quasiparticle 
excitation is of the Bogoliubov type and is given by, 

e k = vVfc) 2 + 2 n U e° (62) 

where e2 = C/n = 47r ^ 2a is the strength of the effective interaction 

between the atoms, n is the density of atoms and m is the mass of the atoms. 
This spectrum is not similar to that of superffuid helium. In superffuid he- 
lium the quasiparticle excitations exhibit roton minimum because of strong 
short range correlation between the helium atoms, but due to the diluteness 
of the atomic vapor quasiparticle excitation spectrum does not exhibit roton 
minimum. The Bogoliubov spectrum is linear in the low wavenumber limit 
and quadratic in the large wavenumber limit. The transition from the linear 
to quadratic spectrum occurs when e° k ~ nU . The reduced wavelength of the 
phonon corresponding to this transition point is called coherence length(£). 
For wavelengths less than the coherent length of the quasiparticle the spec- 
trum is phononic. Since the velocity of sound is c = ^JnUo/m , the coher- 
ent length(£) = h/y2mc. For a singly quantized vortex r ergo = h/mc, so 
27r£ = A = \/2irr ergo . 

(a) BEC of Sodium, 23 Na : In experiments conducted by the MIT group, 
a BEC of sodium atoms was produced with the following parameters: 

• Number density of Sodium atoms, (n) = 4 x 10 20 m -3 . 

• Speed of sound, (c) = 10.4 mm/ s . 

• Coherence length (£) = 0.187 fim. 

Using the above data the following are calculated to be: 

• Wavelength of phonon (X p h) = 1.172 /im. 

• Radius of the ergo sphere, (r ergo ) = 0.263 \xm. 
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• Maximum angle of bending of phonon, S(p = 6.84° 

max 

• Minimum focal length,/ = 4.91 \xm 

min 

• Maximum time delay is, St = 80 /is. 

max 

(b) BEC of Rubidium, 8T Rb : We also have from the JILA experiment us- 
ing Rubidium: 

• Number density of Rubidium atoms, (n) = 2.6 x 10 20 m~ 3 . 

• Speed of sound, (c) = 3.0 mm/s . 

• Coherence length (£) = 0.165 fxm. 

Using this data-set the calculated values are found to be: 

• Wavelength of phonon (\ p h) = 1-037 fxm. 

• Radius of the ergo sphere, (r ergo ) = 0.233 jj,m. 

• Maximum angle of bending of phonon, 50 = 6.84° 

max 

• Minimum focal length,/ = 4.37 \xm 

min 

• Maximum time delay is, St = 0.237 ms. 

max 

Looking at the order of magnitude of the observable, it is evident that BEC 
in atomic gases are more suitable than liquid helium for experimental ver- 
ification of the bending of a phonon trajectory and the corresponding time 
delay. Even though the order of magnitude of the observables for the BEC 
is small, given the recent advances in technology, they are detectable. 
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9 Discussion 



In acoustic analogue models of gravity, propagation of a sound wave in an 
inhomogeneous flow of barotropic inviscid fluid is equivalent to that of a 
minimally coupled massless scalar field in curved acoustic geometry which is 
Lorentzian in signature. In the eikonal limit we can treat a sound wave as 
a ray and also consider only the propagation of a single phonon instead of a 
sound ray in the acoustic geometry. In this limit the trajectory of a phonon 
is a null geodesic in the acoustic geometry. When such a phonon travels 
through an inhomogeneous flow of a barotropic inviscid fluid, its trajectory 
is bent. From the perspective of acoustic geometry this is the bending of the 
trajectory of a massless particle due to the curvature of the acoustic geometry. 
Due to this bending of the trajectory, a phonon has to travel along a path 
which is longer than the path that it would have taken had the fluid been 
at rest. In the short wavelength limit, the leading terms in the expansion 
of the bending angle in powers of the velocity parameters, are quadratic. 
This result is consistent with Fischer and Visser'sjH] result for an irrotational 
vortex without a sink. In the same limit the expansion for the time delay in 
powers of the velocity parameters (A and B) to first order, is independent 
of A. This phenomenon in acoustic geometry is similar to the light bending 
effect and the corresponding time delay in a physical gravitational system 
when light passes by a massive object. 

Instead of a single phonon if we consider two phonons propagating in 
the same direction but on opposite sides of the vortex, they will be seen 
to converge to a point. Thus a vortex in an inviscid, barotropic fluid acts 
like a convergent lens. Both the angle of bending and the time delay in the 
case of a BEC are significantly larger than those for liquid helium. In this 
respect BEC's in atomic gases are more suitable than liquid helium for ex- 
perimental verification of these effects. Since this acoustic analog model of 
gravity is based on the assumption that when a sound wave passes through a 
flowing fluid, the linearized perturbations of dynamical quantities are stable, 
demonstrating the stability of the vortex under such perturbations is more 
important an issue than concerns about the order of the magnitude of the 
observable. This is one of the directions in which this work can be extended 
in future. 
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